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ABSTRACT 

The 3 seismic sources SI, S2 and S3 detected from MDI dopplergrams using 
the time-distance diagram technique are presented with the locations, areas and 
vertical and horizontal velocities of the visible wave displacements. Within the 
datacube of 120 Mm the horizontal velocities and the wave propagation times 
slightly vary from source to source. The momenta and start times measured 
from the TD diagrams in the sources S1-S3 are compared with those delivered to 
the photosphere by different kinds of high energy particles with the parameters 
deduced from hard X-ray and 7-ray emission as well as by the hydrodynamic 
shocks caused by these particles. The energetic protons (power laws combined 
with quasi-thermal ones, or jets) are shown to deliver momentum high enough 
and to form the hydrodynamic shocks deeply in a flaring atmosphere that allows 
them to be delivered to the photosphere through much shorter distances and 
times. Then the seismic waves observed in the sources S2 and S3 can be explained 
by the momenta produced by hydrodynamic shocks which are caused by mixed 
proton beams and jets occurring nearly simultaneously with the third burst of 
hard X-ray (HXR) and 7-ray emission in the loops with footpoints in the locations 
of these sources. The seismic wave in the source SI, delayed by 4 and 2 minutes 
from the first and second HXR bursts, respectively, is likely to be associated 
with a hydrodynamic shock occurring in this loop from precipitation of a very 
powerful and hard electron beam with higher energy cutoff mixed with quasi- 
thermal protons generated by either of these 2 bursts. 

Subject headings: Sumflares — Sun: X-rays, gamma-rays — Sumhydrodynamics 
— Sun: helioseismology 



department of Computing, University of Bradford, Bradford, UK; v.v.zharkova@brad.ac.uk 
2 Department of Applied Mathematics, University of Sheffield, Sheffield, UK; s.zharkov@sheffield.ac.uk 



-2- 



Introduction 



The first successful attempt to observe seismic waves in the SOHO/MDI dopplergrams 
in a form of ripples centred on the X2.6/1B solar flare of 9th July 1996 was reported by 
Kosovichev & Zharkova (1998). The authors presented a helioseismic response (solar quake) 
propagating on the solar surface from the flare location with a strong localised plasma down- 
flow of about 1.5 km/s in the loc ation of the H n flare impulse occurring within a minute 
close to the hard X-ray maximum (jKosovichev fc Zharkova 19981 ). 



A comparison of the observations with the theoretical model (jKosovichev fc Zharkova 19951 ) 
revealed that the momentum required to produce the observed seismic response (~ 2 ■ 10 22 
g cm s _1 ) is one order of magnitude higher th an those of ~ 10 21 q cm s" 1 ob served from the 



plasma downflows in the MDI dopplergrams (jKosovichev &: Zharkova 19981 ). The required 
momentum could be delivered by a hydrodynamic shock appearing at the injection of a very 
hard (7 = 3) and intense (F = 10 12 erg /cm 2 / s) electron beam. These parameters were 
arbitrary selected since there are no detailed hard X-ray observations available for this flare. 

However, the travel time of this shock to the photosphere is more than 2 minutes while 
the time, at which the helioseismic response started in TD diagrams, coincides very closely 
with the time of the hard X-ray impulse and does not reveal the 2 minute delay. Hence, 
there should be some additional sources that can deliver the required momentum to the solar 
photosphere within a very short timescale coinciding with the start time of a hard X-ray 
impulse. 

Recent observations, which reported helioseismic emission from the solar flares of 2003 
October 28 and 29 using the helioseismic holography technique (Donea and Lindsey, 2005; 
DL05 thereafter), revealed another 5 sources occurred in the active region NOAA 10486 with 
the seismic emission at frequencies from 3 mHz to 7 mHz, 4 of them in the flare 28 October 
2003 (DL05) . The flare hard X-ray e mission at the ver y start at 11:01:00 UT was observed by 
KORONAS fcuznetsov et al. 20pd \ and INTEGRAL jGross et al. 20041: iKiener et al. 2006h 
while the RHESSI started observations only at 11:06:00 UT (IHurford et al. 20061 ). However, 
two of the 4 sources for this flare reported by DL05 were wel l aligned with the h ard X-ray 



and 7-ray signatures observed by RHESSI after 11:06:00 UT (IHurford et al. 20061 ) 



Hence, in order to establish a connection between high energy particles and the seismic 
source agents, let us investigate in more detail the velocities of vertical and horizontal dis- 
placements, or the ridges, associated with these seismic waves in the flare 28 October 2003, 
by applying the time-distance diagram technique (TD-method, thereafter) (Kosovichev and 
Zharkova, 1998) to the MDI dopplergrams and by deducing the momenta required to cause 
the observed ridges. Then we can compare them with those delivered by high energy particles 
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of different kinds via low temperature hydrodynamic shocks occurring in the chromosphere 
in response to the particle injections and travelling towards the photosphere. 

The observations used in this study are presented in § [2] describing the active region 
morphology in § 12. 1[ the available high energy observations in § 12. 2\ the time-distance di- 
agrams technique in § 12.31 and the observed ridges of the seismic sources detected in § 12.41 
A theoretical basis for the energy transport from the corona to the lower atmosphere is 
described in § [3] including the heating functions of different particles in § 13.11 and a hydrody- 
namic response to the injection of particle beams or jets in § 13.21 The evaluation of the wave 
parameters in the helioseismic sources and their comparison with the momenta delivered by 
hydrodynamic shocks for different kinds of particles is discussed in § HI The conclusions are 
drawn in § 



2. Description of the observations 
2.1. Active region morphology and magnetic field 



The X17.2 flare occurred on the 28th October 2003 in the very active region NOAA 
10486 at the location 18E20S. It started as observed by GOES from 9:41 UT lasting until 
11:24 UT with the maximum at 11:10 UT in soft X-rays and until 18:00 UT in H a emission. 
The active regions NOAA 10486 had a complex delta-sunspot and produced dramatic flare 
activities in the descending phase of the solar cycle 23 with 3 X-class flares, i.e. an X17 flare 
on 2003 October 28, an X10 flare on 2003 October 29, and an X8.3 flare on 2003 November 



2 and many weaker ones (jYan et al. 20051 ). 



By tracing the changes of the sunspot group simultaneously with the TRACE white- 
light images, the penumbral segments are found to decay rapidly and permanently right after 
each of three X-cl ass solar flares occurred in this region with the neighbouring umbral cores 
becoming darker (jLiu et al 20041 ). These variations are concluded to reflect the changes of 
the photospheric magnetic fields associated with the decaying penumbral areas with some 
parts of them being converted into the umbral field. This implies the emergence of a new 
magnetic flux along the magnetic neutral line and a strong magnetic shear developed in this 
active region that plays an impor tant role as the trigger of the X-class flare on October 28 
JYan et al. 20051 : llshii et al. 20041 ) . 



In the present paper we use the MDI dopplergrams obtained aboard SOHO from 11:00 
UT until 12:00 UT with 1 minute cadence supported by the MDI magnetograms and white 
light images taken from the Solar Feature Catalogues (SFCs) at the times closest to the flare 
setup time flZharkova et al. 20051 ). In Figured] all the sunspots detected in the MDI white 



-4- 



light image at 11:05:33UT with the automated technique (IZharkov et al. 20051 ) taken from 
SFCs are depicted with their umbras and pores (the upper plot) and overplotted onto the 
dopplergrams taken at 11:05:33 UT (from Figure [2]) (the lower plot). It appeared that 3 out 
of the 4 seismic sources X1-X4 reported by DL05 occurred either around the new umbras 
appearing next to the old ones (XI, X3) or in a new magnetic flux of the opposite polarity 
appearing in the existing umbra (X4) (see Figure [2]). 



2.2. Hard X-ray, 7-ray emission and accelerated particles 



The flare 28 October 200 3 started in high energ y emission before 1 1:00 UT that was 



only observed by K ORONAS (IKuznetsov et al. 2000 ) and INTEGRAL dGross et al. 2004 



Kiener et al. 20061) w hile the RHESSI payload did not start the observing until 11:06:00 UT 

(IHurford et al. 20061 ). There is also a very strong CME, an interplanetary shock wave with an 

onset time of 11:01:39 UT an d high energy particles registered at the Earth orbit (IKuznetsov et al. 2006 
Miroshnichenko et al. 20051 ). 



The flare light curves in 7-rays measured by the SONG SOlar Neutrons and Gamma 
rays) instrument aboard KORONAS are plotted in Figured with t he particle energy spect ra 
in Figure [3b (a courtesy of Dr. V.Kurt and the KORONAS team (IKuznetsov et al. 20061 )). 
Similar light curves were observed in the 7-continuum (2.8-3.7 MeV and 7.6-10 MeV) and 
7-lines (2.22 MeV, 4.44 MeV an d 6.13 MeV) by the instruments aboard INTEGRAL satellite 
(Figure 5 in ([Gross et al. 20041 )). 



These light curves revealed the three distinct phases in the flare evolution: a short 
impulsive phase A (under one minute) with a sharp increase of the continuum emission in 
both channels with the photons detected up to 15 MeV, the SPI energy limit, and the two 
longer phases B and C, contained a sharp increase of the 2.8-3.7 MeV continuum emission 
(phase B) and much smoother increase of the line emission in all 3 lines (phase C). 

The start time of the flare impulsive phase was abo ut 11 : 01 : 00 -r 11 : 02 : 00 UT (see 
Figure [3D based on the measurements by INTEGRAL jGross et al. 20041 ) and KORONAS 
(IKuznetsov et al. 2006h . From 11:01:39 UT until 11:05:40 UT the brightest outburst in 
energy was in the range of 0.5 -j- AOMeV. While the protons > 100 MeV were only detected 
at about 11:06:00 UT when also RHESSI started to observe. The KORONAS light curves 
show the 2 rather distinguished peaks in 7-ray emission at 11:02:00 UT and 11:03:00 UT in 
the ranges from 0.5 MeV to 41 MeV (corresponding to the phase A and the start of B in 
the INTEGRAL) and another 2 peaks appearing later at about 11:06:00 UT in the lower 
energy ranges of 0.15-4 MeV and 26 MeV - 100 MeV corresponding to the phase B from 
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INTEGRAL. 



In the first phase the 2 7-continuum peaks are concluded to be produced mainly by a 
bremsstrahlung spectrum generated by el ectrons with energies up to 150 MeV with a small 
7-ray increase in the range of 1.5-7 MeV ([Gross et al. 2004 ; iKuznetsov et al. 20061 ) . In this 
phase the 7 line emission was also observed by the KORONAS that indicates a presence of 
protons with the energies > 30 MeV but not higher than 200 MeV because of the absence 
of photons from the 7r-decay process. The photon energy spectra obtained during the first 
phase (see Figure [3b) are single power laws with the spectral indices about 2 at the lo wer 
energy part (< 70 KeV) and 3.5-4 at higher energies (IKuznetsov et al. 20061 ; Kurt 20061 ). 



The second phase, a delayed one, where the other 2 peaks are observed in 7-rays by 
KORONAS, INTEGRAL and RHESSI, has very noticeable plateaus in the energy spectra 
in the range of 25-100 MeV (see the middle and bottom rows in Figure [3b) first appearing at 
about 11:06:10 UT. Also there are higher energy protons > 200MeV appearing after 11:0 6:00 
UT indicated by a presence of 7r-decay photons (IKuznetsov et al. 20061 ; IShare et al 20041 ). 



The images of th e sources of hard X -ray (200-300keV) and 7-ray emission (2.2MeV) 
obtained by RHESSI (IHurford et al. 20061 ) are marked as Gl and G2 on the active region 
image in Figure [2j The images reveal that at least after 11:06:00 UT when RHESSI started 
to observe there were 2 footpoints with the hard X-ray and 7-ray sources, which have slightly 
different spatial locations with the seismic sources located between them. 

The spectral indices of the proton energy spectra deduced from the ratios of 12 C and 
ie O lines observed i n the phase B, or after 11:06:00 UT, from INTEGRAL vary from 3 to 3.8 
(IKiener et al. 20061 ). This is close to those of 2.8 =F 0.4 reported also after 11:06:00 UT by 
the RHESSI measurements from the de-excitation line 2.22 MeV and positr on annihilation 
line 511 keV with a total number of protons of about 10 33 ([Share et al 20041 ). 



Hence, in this flare there are ample indications from the observations of high energy 
emission about a few events with particles arriving at various times, starting 11:02:00 UT 
until 11:07:00 UT at the footpoints. 



2.3. The time-distance diagrams technique 



Now let us investigate from the MDI dopplergrams this flare evolution at the pho- 
tospheric level and below. We use the one minute cadence dopplergrams for the hour 
starting from 11:00:00 UT. Since the typical oscillatio n frequencies associated with quakes 
( jKosovichev &: Zharkova 19981 : iDonea &: Lindsey 20051 ) are higher than the background os- 
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cillations (3mHz), we also apply frequency filtering centred at 5-6 mHz in order to increase 
the signal-to-noise ratio. 

The obtained velocity distributions in the area with radius of 120 Mm are fit by a 
circular wave using the angular Fourier transform for the a ngle 9 that can be described by 
the angular Fourier transform (jChristensen-Dalsgaard 20031 ) for the angle 9 as follows: 



v(r,9,t) 



=0,2f m l 



r, t)e 



imO 



(1) 



where m = denotes a circular wave, m — 1 - a di polar wave and m = 2 a a qua druple 
wave. Practically, as the previous observation shown (jKosovichev fc Zharkova 19981 ). it is 
not expected to observe the waves higher than these three types because of the observational 
noise. In the present study we extract only a circular wave, while the dipolar wave was also 
registered for this flare by Kosovichev (2006). 

The MDI dopplergrams are re-mapped into polar c oordinates with the cen tres around 
the location of the holographic seismic sources X1-X4 (IDonea fc Lindsey 20051 ). The new 
velocity distributions v{r,9,t) are obtained as the function of time t and their distance r 
from the centre (defined with the accuracy of a single pixel). 

Then we extract the velocities v®(r, t) for every r averaged over various angles 9 for m = 1 
in the datacube of 120 Mm and measure for the different times the horizontal displacements 
of the propagating seismic waves (jChristensen-Dalsgaard 20031 ). i.e. v\(rx,t\), ^0(^2,^2) that 
can be plotted for different tim es (axis Y) and distances (axis X) as a time-distance diagram 
( jKosovichev fc Zharkova 19981 ) . 



2.4. The observed seismic sources 

In order to precisely detect the seismic wave centres, we selected the areas of 20x20 pixels 
around the locations provided by DL05 (Table 1, the last three columns). Then for their 4 
sources (X1-X4) from Table 1 we obtained 400 TD-diagrams, which were visually investigated 
for the pixel locations with detectable ridges denoting the seismic wave propagation. Then 
the pixels with the detectable ridges were used to define the total areas for each quake, and 
their centres of gravity were used as the centres of the seismic waves. 

For the detection of seismic sources we use the one minute cadence dopplergrams start- 
ing from 11:00:00 UT until 13:00:00 UT. In general, in all the MDI dopplergrams we have 
detected 11 locations with downward motion larger than 1 km/s (Figure H]) using the au- 
tomated technique (Zharkov et a!., 2005), while only three of them, SI, S2 and S3, have 
revealed detectable ridges, or quakes (Figure [6]). For each seismic source we extracted the 
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datacube of 120x120 Mm centred in the centre of gravity of the downward Doppler motions 
in this source, re-mapped them into the polar coordinates and applied to the Fourier tech- 
nique (formula 1) for a distance from the centre where the initial impulse is applied (the 
center of gravity) up to 120 Mm. The deduced locations for the seismic sources SI, S2 and 
S3 with the detectable ridges are summarised in the Table 1 (columns 1-3) compared with 
the sources X1-X4 detected by DL05 (columns 4-6). 

The downward velocities V vert (Figure ED directly measured from the dopplergrams in 
these 3 sources were 2.15 km/s for SI, 2.0 km/s for S2 and 1.75 km/s for S3, while for 
another 8 of 11 sources the downward velocities were lower than 1.6 km/s . The durations 
of the downward motions do not exceed 1.5 -2.0 minutes, starting at or after 11:05:00 UT 
with the maximum a minute later and then decreasing for another 30-60 seconds back to the 
pre-flare magnitude. 

The time-distance diagrams for our sources S1-S3 are presented in Figure [6] without and 
with the theoretical ray paths (white solid lines) for the source SI (upper plots), S2 (middle 
plots) and S3 (bottom plots). The start times of the seismic waves are slightly different in 
each source varying from about 11:05:00-11:06:00 UT for the source SI to 11:06:00 UT for 
the sources S2 and S3. These times are more than 3 minutes later than the first maximum 
of hard X-rays for the source SI but close to the second maximum in both HXR and 7- 
rays pointing to the presence of high en ergy protons (> 200 MeV) for the sources S2, S3 



flKuznetsov et al. 20061 : IShare et al 20041 ) 



The source SI (Figure [6], upper plots), corresponding to the source X3 (DL05), at the 
very start of the seismic response, or ridge, reveals the initial horizontal velocity of about 42 
km/s that ncreased in 50 minutes up to 140 km/s at the distances of 120 Mm. The ridge 
in the source S2 (middle plots), or the source XI in DL05, was slightly steeper compared 
to the source SI, i.e. the seismic wave was slower, its velocity varies from 38 to 128 km/s 
in 53 minutes. The ridge in the source S3 (Figure EJ the bottom plots), the source X4 in 
DL05, shows the lowest velocity variations from 34 to 114 km/s in 60 minutes. Hence, the 
initial velocity in the source SI is higher and the seismic wave propagates towards the 120 
Mm edge of the datacube slightly faster than in other two sources S2 and S3. 

The areas of the seismic sources, defined by a presence of the downward motions and 
by detectable ridges, were about 5.05 x 10 17 cm? for the source SI, 3.34 x 10 17 cm? for the 
source S2 and 3.22 x 10 17 cm? for the source S3. Then by using the downward velocities above 
(see Table 2) and comparing with the theoretical seismic ridges (Kosovichev and Zharkova, 
1995, 1998) one can deduce the momenta required to cause the observed ridges: 4.0 x 10 22 
g-cm- s- 1 (SI), 3.7 x 10 22 g ■ cm ■ s' 1 (S2) and 3.1 x 10 22 g ■ cm ■ s^ 1 (S3). 
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The locations of the seismic sources SI, S2 and S3 found from the TD diagrams are 
slightly different than those reported from the holographic method DL05 that can be a 
result of different sensitivity of the techniques applied by us and DL05. The TD technique 
does not produce distinguished ridges for the source X2, which was seen after 11:07:00 UT 
(DL05). There was a weak downward source in the location of the source X2 observed by 
MDI between the times 11:02:00 and 11:03:00 UT (the first maximum in hard X-rays) as 
was spotted in Figure 5 by DL05. While the locations of the sources S1-S3 coincide with 
the dark spots inside the sunspot umbras, there is no indication of any new umbras in the 
location of the source X2. 

The absence of the 4-th source in the TD diagrams can have a number of explanations. 
The source X2 can be rather weak and, thus, cannot be detected by the TD technique while 
seen by the holographic method. Other options are either that this source is an interference 
of the seismic waves produced by the two sources (possibly, SI and S2) or it is located so 
closely to the source SI that we merged them into the extended SI source. However, the 
problem with the latter is that the TD diagram for SI shows its start time at 11:05:00 
UT. These are very puzzling questions that we with the authors of the holographic method 
(DL05) are planning to investigate in the future for this and other flares. 



3. The particle energy transport 

In order to evaluate the transfer into the ambient plasma of the particle momenta and 
energy after their injection from the top into a flaring atmosphere, let us investigate their 
heating functions and resulting hydrodynamic responses. 

Let us assume that the protons or electrons accelerated in a reconnecting current sheet 
(RCS) with a strong longitudinal magnetic field occurred during this flare then they can 
be ejected as power law b eams and completely or partially separated into the opposit e 



footpoints of the same loop (IZharkova &: Gordovskyy 2004 : IZharkova fc Gordovskyy 2005al ). 



In addition to the electron or proton beams, let as consider the particles with Maxwellian 
(thermal) energy distributions shifted to higher energies accelerated along the separatrices 
of the RCS (Gordovskyy et al., 2005). 

Therefore, we consider the 4 kinds of high-energy particles whose energy losses are 
converted into the ambient plasma heating: fast electron beams and fast proton beams with 
power law energy distributions and slow electrons and protons of the separatrix jets with 
thermal energy distributions. 
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3.1. The heating functions by high energy particles 

The energy deposition functions, or heating rates, for these particles are calculated 
using distribution functions found from the full kinetic approach by solving the Fokker- 
Plank equation for electrons losing their energy in collisions and Ohmic heating (Zharkova 
and Gordovskyy, 2005b) and protons in the generation of kinetic Alfven waves (KAWs) and 
their dissipation via Cherenkov's resonance (Gordovskyy, 2005; Gordovskyy et al., 2005). 

The volume heating rates by all kinds of high energy particles simultaneously present 
in the flaring atmosphere are calculated from the particle distribution functions as a vertical 
gradient of the beam energy flux: 



where £ = [ Zmax n(z)dz is a column density, i.e. a number of the ambient particles in 
the area of 1cm 2 on a line of sight from the height z m in to z m ax, n(z) = n(£) is a total 
density of the ambient plasma at a given height, F e (£), F p (£) and Fj(£) are the energy fluxes 
carried by fast electrons, "fast" protons and "slow" protons (of separatrix jets), respectively. 
We do not include the energy losses by slow electrons since these are negligable compared to 
the other particles (Gordovskyy et al., 2005). The heating rate per particle of the ambient 
plasma is related to the volume heating rate S(£) as = ^j|y. The variations of the 
density are considered from a hydrodynamic response below. 



3.2. Hydrodynamic response to the particle injection 

Let us now consider a hydrodynamic response of the ID solar atmosphere to the injection 
of electrons and/or protons by taking into account the continuity, momentum and energy 
equations for the ambient electrons and protons/ions. 

The physical conditions in a flaring atmosphere can be described by a plasma density n; 
electron T e and ion Tj temperatures and a vertical velocity v. All these parameters vary with 
a vertical coordinate z, or a column density £, and time t. The ambient plasma response to 
the injection of high energy particles is described by the hydrodynamic equations (see e.g. 
Somov et al., 1981; Fisher et al., 1985a-c): 



a) Continuity equation: 
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b) Momentum equation: 



<9u 1 9 19 <9i> 

_ + __[„ MJ , + lTe)] = + 9s , (4) 



c) Energy equation for ions: 



ksTi— = -rjiti (— ) + Q(n,T e ,Ti), (5) 



7 - 1 ^ <9t 3 ' V S£ 

d) Energy equation for electrons 



7-1 ft e dt ^<9t dC d£ 

Here 7$ and T e are the ion and electron temperatures, respectively, n is the ambient 
plasma density, x is the ionization degree of the ambient plasma, 7 is the adiabatic constant, 
ks is the Boltzmann constant, k is the thermal conductivity, rji is the ion viscosity, \ is 
the full ionization energy of a hydrogen atom, /1 = 1.44m#, g is the acceleration due to 
gravity of the Sun. P(n,£) indicate the volume heating rates provided by electrons owing to 
collisions, by electrons owing to Ohmic losses, by protons owing to collisions and by KAWs, 
L ra d is the volume radiative energy loss rate and, finally, Q(n,T e ,Ti) is the rate of energy 
exchange between ions and electrons. 

The solution is sought in a limited region of the solar atmosphere £ m j n < £ < ^ max with 
the minimum boundary located at £ m j n = 2 x 10 17 cm~ 2 and the maximum boundary is 
located deep in the photosphere at £ max = 2 x 10 24 cm~ 2 . The initial atmosphere is assumed 
to be in hydrostatic equilibrium, i.e. v(0, £) = and isothermal, i.e. T e (0, £) = Tj(0, £) = T 
where T = 6700° K. The ionization degree x is defined by a modified Saha formula (Somov 
et al., 19 81). We also take int o account the initial momenta delivered by the particles at 



injection (IBrown fc Craig 1984 ) by assuming v(0, £ m i„) = — where Pb is the sum of the 



momenta of all injected particles. 

The initial distribution of a plasma density is defined as follows: 



n(0,f) = n min + Iiq 1 ^ - £ min ), (7) 
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where n m i n = 10 cm and ho is the height scale: 

A = + (8) 

Me 

The radiative losses rate L(n,T e ) is described by the analytical expression: 



L n ,T e = n 2 xL(T e ) + nL H (n,T), (erg/cm 3 /s) (9) 
where the radiative l oss function L(T e ) is taken for the coronal abundances of elements 



in optically thin plasma (ICox fc Tucker 19711 ) and Lfj(n, T e ) are the radiative losses in a ll 



hydrogen lines calculated for the optically thick atmosphere (IZharkova fc Kobylinskij 19931 ) . 



The boundary conditions are defined as follows. 1). We assumed that there is no an 
initial heat flux on the top boundary, i.e. dTe( -°^ min '> = gZM°i|™rb> = 0; 2) the upper boundary 
is a free surface in the presence of the coronal pressure, i.e. 

dv(t,j; min ) = 4J_^ ^ -Pcor(z(t,£min))], (10) 

where p(t,£ min ) = nk B {Ti + xT e ) and p C or{z{0,t, min )) = n min k B [l + x(T Q )]T where the 
ionization degree x is defined by a modified Saha formula (Somov et al., 1981). 



3.3. The momenta delivered by beams and hydrodynamic shocks 

3.3.1. The momentum delivered by a proton and/ or electron beam 
The momentum P SjP delivered in pure collisions by an electron o r proton beam with a 



spectral index 7 e iP and a lower energy cutoff Ei OWe can be evaluated as (IZharkova & Gordovskyy 2005bl ) 

/•OO 

P e , p = ^/2m~K / E~^ +0 - 5 dE, (11) 

J Elow e ,p 

where m e p is the electron or proton mass, E\ is the lowest energy in the relevant particle 
spectrum (theoretical lower energy cutoff) and K is the normalisation constant that can be 
found from the total number of measured particles N e ^ p : 



OO TV- 

n 7 + 1 



N e , P = K / E'^dE = -E-i; v \ (12) 
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where Ei OWep is the measured lowest energy of electrons or protons. Hence, without 
taking into account pitch angle scattering, the momentum delivered by a proton beam can 
be evaluated by substituting the constant K found from N p into the equation for P ep and 
performing the integration. This will result in the following: 

£-(7e,„-0.5) 

P e>p ~ y^T p N e , p gg- 1 — . (13) 

( 7e , p + 0.5) E j£* D 

If we assume that E\ = Ei OWep , then the momentum delivered by electrons/protons without 
pitch angle scattering can be evaluated as: 



7. 



P e , p ~ ^2m^ p N e , p ■ Ej™. (14) 



7e,p + 0.5 



It should be emphasized that this is the upper limit of the momentum carried downwards 
to the photosphere by electrons or protons since it is calculated without taking into account 
pitch angle scattering and wave dissipation fo r protons or Ohmic diss ipation for electrons 



that can reduce its magnitude by a few factors (iGordovskyy et al. 20051 ). However, it allows 
the comparison of the momenta delivered by high energy particles with those measured from 
the downward motions in dopplergrams and from the TD diagrams. 



3.3.2. The momentum delivered by a hydrodynamic shock 

Let us also evaluate the momentum delivered by a hydrodynamic shock using the simple 
formula: 

P hd = E t mv(t) (15) 

where the summation is done over the time from to r where r is a duration of the 
impact causing the seismic waves, m is the mass of the plasma delivering the momentum 
related to the flaring area A where the momentum is deposited, V is a starting velocity of 
the plasma at the moment of impact and t is a duration of the impact. 

For the known plasma mass density p = m# • n where n is the particle density per 
volume defined from hydrodynamic solutions, this equation can be re-written as follows: 



P hd = Z t mv{t) w p-A-v 2 • r, 



(16) 



13 



where p is an average density of the plasma delivering the momentum, A is the flaring 
area where the momentum is deposited, v is an averaged velocity of the plasma propagation 
at the moment of impact and r is the duration time of the impact. 



4. Results and discussion 

4.1. The momenta delivered by beams 

Let us try to establish the agents delivering the momenta reported in Table 2 (§ 12.41) by 
investigating the parameters of high energy particles associated with each source from the 
hard X-ray and 7-ray observations by CORONAS, INTEGRAL and RHESSI in § E2J 

The source SI 

For the source SI one has no 7-ray emission but only hard X-ray photon differential 
spectra (or mean flux) as presented in Figure [3b with the upper energy of 150 MeV and the 
lower energy cutoff about 18 keV (the left upper plot for 11 : 02 : 00-=- 11 : 03 : 00 UT, Kurt, 
2006). The photon spectral index in this plot was about Sfugh — 3.5 — 4.0 at the energy 
range from 70 keV to 60 MeV and 5i ow = 1.5 for the energy range of 10 -r 70 keV (Figure [3]). 
Since this was a very strong flare, we assume that the precipitating beam has induced a very 
strong; electric field and its hard X-ray emission was dominated by the Ohmic energy losses 



( IZharkova fc Gordovskyy 2005bl ; IZharkova fc Gordovskyy 20061 ). This assumes the spectral 
index 7 of a precipitating electron beam has to be nearly the same as the index 5hi g h of the 
photon spectrum at higher energy. 

The difference between the spectral indices Shigh an d 6i ow can provide us with the beam 
initial energy flux F for the selected spectral index of an electron beam, i.e. for the beams 
with 7 varying from 3. 5 to 5, the difference in the ph oton indecies varies from 2.0 to 2.5. 



Hence, from Figure 12 (IZharkova &: Gordovskyy 2005bl ) we can deduce the initial energy flux 



of bea m electrons that can vary from 1 ■ 10 erg /cm / s for 7 = 3.5 to 4 ■ 10 erg /cm / s for 



7 = 5 (IZharkova &: Gordovskyy 2005bl ). These fluxes can be carried out by electron beams 



with the initial densities of 1 • 10 cm and 6 ■ 10 cm , respectively. 

Let us calculate the momenta delivered to the photosphere by electron beams with such 
parameters using the technique described in § 13.3.11 For the flare area of about 5.05 x 10 17 cm? 
defined by the area of a downward Doppler motion Figure H] in the source SI, the momentum 
delivered to this area by the electron beam was about 2 • 10 20 g ■ cm/s. Obviously, this is 
not sufficient to deliver the required momentum of (3 — 4) x 10 22 as reported in Table 2 and 
to directly cause the ridge observed in the TD diagram for SI (see § 12. 4p . 
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The sources S2 and S3. 

The sources S2 and S3 appear close to the locations o f hard X-ray emissio n and within 



the circle denoting 7-ray emission observed by RHESSI (IHurford et al. 20061 ). As it was 
noted in § 12.21 the spectral indices of the proton energy spectra observed by INTEGRAL 
in t he phase B after 1 1:06:00 UT (from the ratios of l2 C and 16 lines) vary from 3 to 



3.8 (IKiener et al. 20061 ) or by the RHESSI measurements (from the de-excitation line 2.22 



MeV and positron annihilation line 511 keV) vary 2 .8 =F 0.4 with a total number of protons 



observed estimated at about 10 33 ( jShare et al 20041 ) and the lowest energy about 30MeV, 
or pd 4.8 • 10~ 5 erg. However, this energy can be decreased to 2 MeV without affecting the 



2.22 MeV and annihilation line emission (jShare et al 20041 ) . 



Then the momentum delivered by such a proton beam to the chromosphere where the 
7-emission is measured, can be about ~ 2.2 • 10 22 g ■ cm/s for the lower energy of 30 MeV 
and according to equations (12) and (14), increases as (2MeV/30MeV)^' r+0 ' 5 x ^ Q 5 to 
~ 5.2 • 10 24 g ■ cm/s for the lower energy of 2.0 MeV and 7 = 3.0. This range superbly covers 
those momenta derived in the sources S2 and S3 (Table 2). Hence, protons can be the agents 
in these two sources delivering the sufficient momenta to the region where MDI measures the 
Doppler velocities. Of cource, as pointed in § 13.3.11 the momentum evaluated from formulae 
(12-14) provide the upper limit since other proton or electrons scattering mechanisms can 
slightly reduce a number of particles reaching a given depth from the top depending on beam 
parameters. However, for moderate hard beams, as reported for this flare (spectral index of 
3.3), this difference is not very noticeable (Zharkova and Gordovskyy, 2005b). 



4.2. Simulated heating functions 

Now let us compare the heating rates of the three kinds of particles considered in § [3] 
and § 13 . 11 The timescale within which particle of each kind can reach the photosphere is 
about 1 s for electron beams, 2-5 s for proton beams and 10-20 s for protons of separatrix 
jets for the standard loop length of 10 9 cm (Zharkova and Gordovskyy, 2004; Gordovskyy 
et al. 2005). Therefore, the propagation time for each kind of particles is short enough 
to contribute to the ambient plasma heating and to form a shock, or a lower-temperature 
condensation, appearing as a result of the hydrodynamic response to these particle injection. 

The heating rates simulated from the full kinetic approach are presented in Figure [7J 
curve A - for an electron beam with the initial energy flux Fq = 1.4 ■ 10 8 erg/cm 2 /s, spectral 
index 7 = 2 and the lower energy cut-off lQkeV, curve B - for the separatrix jet protons 
with the initial energies E = lMeV and initial energy flux of 4 • 10 u erg/cm 2 /s , curve C - 
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for an electron beam with the initial energy flux F = 10 w erg /cm 2 / s, spectral index 7 = 5 
and the lower energy cut-off 16 keV, curve D - for a proton beam with the initial energy flux 
F Q = 4- 10 10 erg/cm 2 /s, spectral index 7 = 1.5 and the lower energy cut-off 40 keV and curve 
B for the separatrix jet protons with the initial energies E = lMeV and initial energy flux 
of 4 • 10 11 erg/cm 2 /s. 

It can be seen that the heating by electron or proton beams with power law energy 
distributions is strongly dependent on the initial beam parameters: softer and weaker beams 
deposit their energy mainly in the corona and upper chromosphere while harder and more 
powerful beams deposit more energy deeper in the lower chromosphere (compare the curves 
A and C for hard and soft electron beams and D for a soft proton beam) (Zharkova an 
Gordovskyy, 2005b, Gordovskyy et al, 2005). 

Electron beams are considered to deposit their energy more evenly in depth compared 
to the proton ones (compare the curves A,C with D). Proton beams deposit the bulk of their 
energy via generation of KAWs with their following dissipation in Cherenkov's resonance 
at the flaring atmosphere depths where their velocities are higher than the local Alfven 
ones (Gordovskyy et al., 2005). It can be noted that the heating by KAWs induced by 
protons has the two regions where this condition stands: in the upper corona because of 
their initial exponential distributions (the first curve B) and at the lower chromosphere 
because of reduction of the local Alfven speed and of the proton exponential distributions 
(the second, spike-like curve B). 

The heating of the upper atmosphere by proton beams can be even more noticeable after 
Coulomb collisional losses are taken into account which will make the proton distributions 
in the chromosphere even more exponential. While strongly affecting the heating by proton 
or electron beams of the corona and upper chromosphere (before the column depths of 
10 20 cm -2 ), the collisions do not change significantly the heating of the lower chromosphere 
where the Cherenkov resonance is dominant (Gordovskyy et al., 2005). Another heating 
mechanism considered for beam electrons is Ohmic dissipation in a self-induced electric field 
that contributes significantly to the heating of the coronal levels but again has not affected 
the lower atmosphere heating (Zharkova and Gordovskyy, 2005b). 

The effect of these heating functions on the hydrodynamic solutions is discussed below. 
We include into the heating functions: collisional losses by both electron and proton beams, 
Ohmic losses only for an electron beam and Cherenkov resonance for thermal-like protons. 
Two heating functions are considered: by a pure electron beam and by a proton beam 
combined with quasi-thermal jet protons. 
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4.3. Simulated hydrodynamic responses 

In order to maximise a deposition at lower atmospheric levels let us compare the hydro- 
dynamic responses caused by a hard electron beam with 7 = 3.5 and the intial energy flux of 
10 11 erg /cm 2 /s (that is higher and harder than measured) with those by proton beams/jets 
with the spectral indices of about 3 and the initial energy fluxes of 10 12 erg /cm 2 / 's as deduced 
from the hard X-rays and 7-rays in § 12. 21 

We simulate the hydrodynamics of a two-temperature plasma heated by either electrons 
or protons (power laws plus those with Maxwellian energy distributions) by solving the 
two energy equations (for electrons and protons), continuity and momentum conservation 
equations and include the radiative cooling as described in § 13.21 

The variations of temperature, density and macro-velocity simulated for the hydrody- 
namic responses are plotted for the electron beam with the parameters derived from hard 
X-rays (Figure El the left plots) or for the mixed proton/jet beam with the parameters 
derived from the 7-ray emission (Figure [HJ the right plots). 

The presented hydrodynamic results caused by beam electrons agree with those by So- 
mov et al. (1981), Nagai and Emslie (1984) and Fisher et al. (1985a-c). For the intense hard 
beam as deduced from the hard X-ray data from KORONAS (Kuznetsov et al., 2006) there 
is a noticeable decrease (rarification) of the total plasma density and a strong temperature 
increase in the corona accompanied by explosive evaporation of the chromospheric plasma 
into the corona (macromotion upwards) occurring in response to the beam injection. 

Starting from the column density 2 x 10 19 cm 2 , a collisional stopping depth for lower 
energy electrons of 12 keV accepted in our simulations, a low temperature condensation is 
formed moving as a shock with velocities of about 100-200 km/s, which are higher than the 
local sound velocity (Somov et al., 1981, Fisher et al., 1985). However, such a shock produced 
even by a powerful electron beam with the spectral index of about 3 (the lower plots) appears 
rather high in the upper chromosphere between the column depths of 2 x 10 19 -j- 2 x 10 20 
cm -2 as can be seen in Figured (the left plots). 

These two motions of the ambient plasma (upwards and downwards) are reported in all 
the hydrodynamic simulations (Somov et al., 1981, Nagai and Emslie, 1984, Fisher et al., 
1985 a-c) and widely investigated from the blue and red-shifted spectral measurements in 
UV and H a emission, respectively. For some events, or some beam parameters, these can be 
nearly equal (e.g. see Zarro et al., 1988) while for many other beams if the electron Ohmic 
losses are taken into account (Zharkova and Gordovskyy, 2006) only the blue shifts could be 
observed without the noticeable red ones. 
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As can be seen from Figure [HI the hydrodynamic response to the injection of the pure 
electron or mixed proton/jet beams are substantially different. The electron beam injected 
during 10 s produces a smaller (by factor 2-2.5) temperature increase in the corona (top left 
plots), a smaller (by an order of magnitude) density depression of the coronal plasma into 
the chromosphere (the middle left plots) and smaller (by factor 2-2.5) evaporation velocities 
(the bottom left plots) compared to those by mixed proton/jet beams (the right plots in 
Figure [8] and Figure [9]) . 

In addition, the mixed proton/jet beam forms the lower temperature shock, which is 
spread much deeper into the lower chromosphere between the column depths of 2 x 10 20 -j- 
8 x 10 21 cm -2 (see Figure [9j:or a closer view of the shock in the first 100 s). The velocities 
of the shock induced by the mixed proton beam are also higher (by a factor of 2-2.5) than 
those induced by the pure electron beam. These macro-velocities induced by the proton 
beam decrease in the region with a column density of 5 x 10 21 cm~ 2 to a few km/s compared 
to those less than 1 km/s for the one induced by electrons. The momentum induced by 
this shock is transferred to the photosphere within much shorter timescale under 1 minute 
because it is formed in much deeper and denser atmospheric levels. This also is confirmed by 
the temporal profile of the macro- velocity variations at the lower edge of the shock (the very 
right points in the distributions in Figure [9]) induced by proton beams shows an increase, to 
a few km/s, of the edge macrovelocities, within a minute, that resembles those measured in 
Figure 

The proton induced shock deposits its momentum from the depths 5 — 8 x 10 21 cmT 2 to 
very dense plasma beneath that is delievered with a velocity of about 2 km/s through about 
120 km of the solar atmosphere before approaching the column depth of 4 x 10 23 cm~ 2 for 
Ni line region. While the electron-induced shock is required to travel from the column depth 
of (2 — 5) x 10 20 cm~ 2 , or about 350 km, with the same or the twice lower velocity that will 
allow the momentum to reach the Ni region with a delay of 3-6 minutes. 

4.4. The region for Ni 6768 line formation 

In order to understand the effect of the hydrodynamic shocks on the doppler measure- 
ments in Ni 6768A line by the MDI instrument, one needs to establish the line formation 
region. This can be done by using the full non-LTE simulations (without magnetic field 
effects) of the Ni line in the ambient plasma with the coronal abundance of elements and 
molecules (up to 23 in total are considered, the main ones are CO, C2, CH, CN)(Bruls, 
1993; Uitenbroek, 2001; Zharkova and Kosovichev, 2002). The ambient plasma tempera- 
ture, density and macrovelocities are described by the hydrodynamic solutions presented in 
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§£31 

We consider the two cases: the quiet Sun atmosphere outside sunspots (Figure [LTJ 
upper plot) before the beam onset and the flaring atmosphere heated by beam electrons 
(Figure [101 bottom plots). The contribution functions for the regions, where the Ni line 
6768A is formed, are plotted by the grey curves with the grey arrows marking the maximum 
contribution. The background contribution functions (for the continuum) are plotted by the 
blue curves with the blue arrows marking a region of the maximum continuum contribution 
around the Ni line. The total contribution functions for all elements are plotted by the black 
curves. The red curves show the Plank contribution functions and the yellow lines present 
the mean intensity J from a given atmospheric level for the temperature distributions taken 
from Figure [H^the top left plot, time = 7s after the beam injection). 

The formation region for the Ni line 6768A lies approximately within the column depths 
of (2.0 — 4.0) • 10 23 cm -2 , or around 200 km for the quiet Sun (Figure [TUl upper plot) that 
increases to (4.0 — 6.0) • 10 23 cm -2 , or around 180 km, for the flaring atmosphere heated by 
beam electrons (Figure [JTJ bottom plot) that agrees rather well with the estimations found 
from the previous non-LTE simulations (Jones, 1989; Bruls, 1993; Zharkova and Kosovichev, 
2002). 

Let us compare these column depths with those obtained for the lower temperature 
shocks simulated from the hydrodynamic responses to injection of the electron or pro- 
ton/electron beams as discussed in § 14.31 One can notice that the column depths for the Ni 
line region are closer to a hydrodynamic shock formed by the mixed proton/jet beam (Figure 
[9] and Figure [HI right plots). Also the shock velocities decrease to a few km/s towards the 
deeper depth edge (the right ends of each curve in Figured]) and this decrease has a temporal 
profile of the lower edge velocities increasing from 0.1 km/s (at Is) to 2 km/s (at 50s) and 
then decreasing back similarly to those measured by MDI (see Figure ED- 

As discussed above in § 14.31 this shock needs only about 60 s to reach the Ni region 
that allows the detection of the seismic response nearly simultaneously with the hard X-rays 
as reported for this flare in § [2] and § 12.41 While the shock caused by the electron beam is 
formed much higher in a flaring atmosphere (Figure El left plots), it can be measured in the 
Ni 6768A line delayed by up to 6 min compared to the time of X-ray and 7-ray emission. 

4.5. The momenta delivered by shocks and the proposed scenarios 

In § !3.2l we established that the hydrodynamic shock, or lower temperature condensation 
induced by the proton beams is denser and wider in depth moving twice as fast as those 
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induced by the electron beam. The travel time to the photosphere of the shock caused by 
electron beams is about 180-360 s (compare the depth and densities of the lower-temperature 
condensations plotted in the lower graphs, in Figure [8]) before it reaches the photosphere. 

The hard proton beam of a moderate power (10 12 erg /cm 2 /s mixed with the thermal 
jet protons presented in Figure [8] and Figure [9]) forms a shock much deeper in the lower 
chromosphere than the electron beam and it travels to the photosphere in about 60 s. The 
shorter travel time of such a shock can explain the close correlation of the onset times for 
the seismic waves, measured in the TD diagrams (Figure [6]) with the hard X-rays and 7-ray 
emission bursts reported in § 12.41 

In order to confirm it, let us calculate the momenta delivered by the HD-shocks caused 
by either beam by using the method described in § 13.3.21 and taking the particle densities, 
macro-velocities and duration times from Figure [Hand the source areas from Table 2. 

The shock caused by a pure electron beam (Figure [HJ left plots) with the beam param- 
eters taken from Figure [3] has a density about 5 x 10 12 cm -3 , an average macro-velocity 
about 1.8 x 10 7 cm/s and a duration about 100s. This shock is formed at the depth of about 
(2 — 5) x 10 20 cmT 2 . The momentum still requires a timescale > 180s to be delivered to 
the Ni line formation region (2 — 6) x 10 23 cm -2 . Then for the source SI with the area of 
about 5.05 x 10 17 cm 2 the electron-formed shock can deliver a momentum of about 1.3 • 10 19 
g ■ cm/s. 

Evidently, the shock produced by a pure electron beam does not contain enough mo- 
mentum to account for the seismic responses recorded in either source. Also a time delay 
of about 180-360 s is required for this shock to reach the region of the Ni line formation to 
cause a delay in the seismic response appearance compared to the emission in hard X-rays 
and 7-rays. 

On the contrary, the parameters of the proton-formed shock are much more relevant to 
the explanation of the momenta observed in the seismic sources SI, S2 and S3. This shock 
occurs much closer to the region where the Ni line is formed, it contains much denser material 
(up to 10 14 cm~ 3 ) and higher macro- velocities ((2 — 3) x 10 7 cm/s. By substituting these 
parameters and the areas of the sources SI, S2 and S3 taken from Table 2 into the formula 
(18), one can find that this shock can deliver the much higher(> 4.4 x 10 22 gem/ s) momentum 
compared to that delivered by the electron-formed one. The momentum induced by this 
shock is transferred to the photosphere within much shorter timescale under 1 minute because 
it is formed in much deeper and denser atmospheric levels. Therefore, the magnitudes of the 
momentum carried by the shock caused by the mixed proton/jet beam and the timescale, 
within which it reaches the region of Ni line, are rather close to those deduced from the 
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TD-diagrams discussed in § 14.11 (Table 2) for all the seismic sources. 

Since 7-ray emission by high energy protons was observed by KORONAS only after 
11:06:00 UT (Kuznetsov et al., 2006), when the sources S2 and S3 appeared, one can assume 
that these sources are produced by a high energy proton beam combined with the lower 
energy (< 200Me^) quasi-thermal protons of the separatrix jets. They together can deliver 
to the photosphere the momenta via the hydrodynamic shock, to cause directly the observed 
seismic responses and to energise ambient electrons to energies high enough for the hard 
X-ray and 7-ray emission observed at 11:06:00 UT in the third burst as discussed in § 12.21 

For the source SI sufficient momentum can be delivered by a very high energy electron 
beam, which is reported to have energies of hundreds of MeVs (Kuznetsov et al., 2006) 
and occurred at 11:02:00 UT in the first burst (Figure 3) combined with the lower energy 
(< 200Mel^) quasi-thermal jet protons. Then the hydrodynamic response to such heating 
could lead to a shock formed slightly deeper in the flaring atmosphere than those for the 
electron beam from Figure [8] and slightly higher than for the protons. Within 3-4 minutes 
after the beam onset the momentum caused by this hydrodynamic shock can reach the Ni 
region and cause the observed seismic waves as discussed above which explains the delay of 
about 3-4 minutes between the seismic response in the source 1 and the hard X-ray emission 
in the first burst started at 11:02:00 UT. 



5. Conclusions 

In the current paper we report the 11 sources with downward motions higher than 1 
km/s in the flare 28 October 2003 and the seismic waves in 3 of them (S1-S3) detected with 
the time-distance (TD) diagram technique. The 3 seismic sources started around 11:05:00- 
11:06:00 UT, had slightly different downward (vertical) velocities and heights of ridges, or 
horizontal velocities, pointing out to different agents causing them. 

We investigate a few agents able to deliver the required energy and momentum to the 
flaring atmosphere: power law electron beam, power law proton beam and quasi-thermal 
protons with energies below 200 MeV occurring through acceleration by a super-Dreicer 
electric field in reconnecting current sheets on the top of a flaring atmosphere (Zharkova and 
Gordovskyy, 2004, 2005a). Electron beams are considered to deposit their energy in collisions 
and Ohmic dissipation, proton beams - in collisions and generation of kinetic Alfven waves 
with their following dissipation in Cherenkov's resonance (Gordovskyy et al., 2005). 

The hydrodynamic responses to a precipitating power law electron beam or to a proton 
beam mixed with the protons of separatrix jets are also investigated. The hydrodynamic 
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response to the injection of the pure electron or mixed proton/jet beams are found to be 
substantially different. The electron beam injected during 10 s produces a smaller (by 
factor 2-2.5) temperature increase in the corona (top left plots), a smaller (by an order of 
magnitude) density depression of the coronal plasma into the chromosphere (the middle left 
plots) and smaller (by factor 2-2.5) evaporation velocities (the bottom left plots) compared 
to those inducedcby mixed proton/jet beams (the right plots in Figure [8] and Figure [9]). 

In addition, the mixed proton/jet beam forms a low temperature shock, which is spread 
much deeper into the lower chromosphere between the column depths of 2 x 10 20 t8x 10 21 
cm~ 2 (see Figure[9]for a closer view of the shock in the first 100 s). The velocities of the shock 
induced by a mixed proton beam are also higher (by a factor of 2-2.5) than those induced 
by a pure electron beam. These macro-velocities induced by the proton beam decrease in 
the region of a column density of 5 x 10 21 cmT 2 to a few km/s compared to those less than 
1 km/s for the one induced by electrons. Also the temporal profile of the macro- velocity 
variations at the lower edge of the shock (the very right points in the distributions in Figure 
|9]) induced by proton beams shows an increase to a few km/s of the edge macrovelocities 
within a minute that resembles those measured in Figure 

The momenta and start times measured from the TD diagrams in the sources S1-S3 
are compared with those delivered to the photosphere by different kinds of high energy 
particles with the parameters deduced from hard X-ray and 7-ray emission as well as by the 
hydrodynamic shocks caused by these particles. 

The energetic protons (power laws combined with Maxwellian ones from the separatrix 
jets) are shown to deliver momenta high enough and to form hydrodynamic shocks much 
deeper in a flaring atmosphere compared to a pure electron beam. This allows the proton- 
formed shocks to travel to the photosphere shorter distances and less time resulting in seismic 
waves occurring nearly simultaneously with the high energy emission as observed in the 
sources S2 and S3. The source SI is likely to be associated with a hard power law electron 
beam mixed with the quasi-thermal protons of the separatrix jets. 
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Table 1: The heliographic locations L and B of the 3 seismic sources detected by the TD 
technique (columns 1-3) and 4 seismic sources i n the flare 28 October 2003 reported from 
the holographic method (IDonea fc Lindsey 20051 ) (columns 4-6). 



Our sources 


L ° 


B ° 


DL05 


L ° 


B ° 


SI 


287.28 


-15.96 


X3 


287.05 


-15.78 


S2 


284.72 


-17.62 


XI 


285.45 


-17.61 


S3 


291.00 


-16.64 


XA 


291.46 


-16.43 








X2 


285.86 


-16.01 
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Table 2: The areas A, downwards (vertical) velocities V vert , the duration of an observed 
downflow motions in the source location, T impu i se , the average horizontal velocity Vhoriz and 
the Momenta deduced for the flare 28 October 2003 from the time-distance diagram for the 
3 seismic sources SI, S2, S3 (Figure [6]). 



Our sources 


area, A, cm 2 


V vert , Cm Is , 


T- impulse-, S 


Vhoriz, cm/s 


Momenta, g ■ cm/s 


SI 


5.05 x 10 17 


2.15 x 10 5 


90 


40.8 x 10 5 


4.1 x 10 22 


S2 


3.34 x 10 17 


2.00 x 10 5 


70 


38.0 x 10 5 


3.8 x 10 22 


S3 


2.22 x 10 17 


1.75 x 10 5 


70 


35.4 x 10 5 


3.1 x 10 22 
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280 285 290 295 
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Fig. 1. — The sunspot group located in the active region NOAA 10486 with the large trailing 
and 2 smaller leading sunspots obtained in the white light image at 11:05:11 UT from the 
Solar Feature Catalogue. 



-28- 



SOHO/MDI 28-Oct-2003 11:1 1:33 



-8.240 




10.240 




12.240 




14.240 




1 6.240 


: J* 


18.240 




20.240 




22.240 
24.240 




26.240 




28.240 






■a - 



1500.0 



1000.0 



500.0 



275 



760 279.760 283.760 287.760 291.760 295.760 
Carr Longitude 

SOHO/MDI 2003-1 0-28T1 1 :05:33 



-8.240 


_ i 


10.240 




12.240 




14.240 




1 6.240 


■a* 




D I -1 


18.240 




20.240 


; V 


22.240 








24.240 




26.240 




28.240 






■a - 



-500.0 



■1000.0 



500.0 



275.760 279.760 283.760 287.760 291.760 295.760 
Carr Longitude 



Fig. 2. — The same sunspot group as in Figure 1 over plotted onto the MDI magnetogram 
11:11:33 UT (upper plot) and on the dopplergram (bottom plot) with the locations of the 
seismic s ources marked by X1-X 4 (columns 4-6 in Table 1) detected by the holographic 
method (IDonea fc Lindsev 20051) w ith the gamma-ray sources observed by RHESSI marked 
as Gl and G2 jHurford et al. 20o3 ). 
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Fig. 3. — The light curves (the two left columns) and the differential spectra (the righ two 
columns) in different energy bands obtained by the SONG (Solar Neutron and Gamma-ray) 
instrument aboard KORONAS during the whol e flare duration (a cou rtesy of Dr. V.Kurt 
(Moscow University) and the KORONAS team (IKuznetsov et al. 20061 ) ). 



-30- 



-7.1 



-10.3 



SOHO/MDI dopplergram 



1000 



'3.5 



I9.8 



m 



-23.0 

220.9 224.1 227.2 230.4 233.6 236. 

Carr. Longitude 

SOHO/MDI dopplergram 11:06 2003-10-28 



-8 


24 


10 


24 


12 


24 


1 4 


24 


1 6 


24 


18 


24 


20 


2-1 


22 


24 


24 


24 


26 


24 


28 


24 



S3 



S2 



- i i 

275.76 278.76 281.76 284.76 287.76 290.76 293. It 
Carrington Longitude 



-1000 



-2000 



-3000 



Fig. 4. — The locations of 11 Doppler sources (red contours) with the downward motion 
higher than 1 km/s (top image) but only 3 of them with the detectable seismic responses in 
the locations specified in columns 1-3 of Table 1 (bottom image). 



Fig. 5. — The downward velocities (m/s) (Y-axis) versus time (minutes) (X-axis) in SI 
source (upper right plot), S2 (upper left plot) and S3 (bottom left plot) in the centre of 
gravity locations with the downward Doppler motions depicted in Figure Hb and Table 1. 
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Fig. 6. — The frequency-filtered time-distance diagrams for the seismic sources: SI without 
a ray path (left upper plot) and SI with the ray path (right upper plot ), similar for the 
sources S2 (middle plots) and S3 (bottom plots). 
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10° 10 1 10 2 10 3 10 4 10 5 
Column density, 10 18 cm -2 

Fig. 7. — The heating rates for different populations of accelerated particles: the curve 'A' 
corresponds to a weak and hard electron beam, the curve "B" marks the heating by "slow" 
jet protons; "C" to a strong and soft electron beam and "D" to a "fast" proton beam. 
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Fig. 8. — The hydro dynamic responses of a flaring atmosphere caused by pure electrons (the 
left column) or by mixed proton and electron beams (the right column) with the parameters 
of protons deduced from 7-rays and electrons from X-rays from RHESSI and KORONAS. 
The upper plots present the ambient electron temperatures, the middle ones - densities and 
the lower plots - macro- velocities, the numbers on the graphs show the times in seconds after 
the beam injection. 
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Fig. 9. — The close-up from Figure [8] of the temporal variations of macro- velocity in cm/s 
in a lower temperature condensation, or a hydrodynamic shock (Y-axis) versus the column 
depth in cm~ 2 (X-axis), appeared in response to the injection of mixed proton/electron 
beams (Figure [HI right plots(l- after 10s (black line), 2 - 30s (grey line), 3 - 50s (purple line), 
4 - 70s (blue line) and 5 - 100s (green line). 
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Fig. 10. — The source function distributions calculated for the quiet atmosphere (top plots) 
and the flaring atmosphere (bottom plots) heated by a hard intense electron beam (7=3, 
F = 10 12 erg /cm 2 / s for the Ni line transition 6768A (S act i ve , grey line), background elements 
{Sbackgr, blue line), Plank function (S p i an k, red dashed line) and total for all elements (S to tai, 
black line) as well as the mean intensity J (yellow line) simulated using the full NLTE 
MULTI-based) approach for the full coronal abundance of elements and some molecules 
(CO, C 2 , CH, CN, 23 in total) (Uitenbroek, 2001; Zharkova and Kosovichev, 2002). 



